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ABSTRACT 


Three models of a flat universe of coupled matter and dark energies with 
different low-redshift parameterizations of the dark energy equation of state are 
considered. The dark energy is assumed to vary with time like the trace of the 
energy-momentum tensor of cosmic matter. In the radiation-dominated era the 
models reduce to standard cosmology. In the matter-dominated era they are, for 
modern values of the cosmologic al parameters, consi stent with data from SNe 
la searches and with the data of iGurvits et all fjl999l) for angular sizes of ultra 
compact radio sources. We hnd that the angular size-redshift tests for our models 
offer a higher statistical conhdence than that based on SNe la data. A comparison 
of our results with a recent revised analysis of angular size-redshift legacy data 
is made,and the implications of our models with optimized relativistic beaming 
in the radio sources is discussed. In particular we hnd that relativistic beaming 
implies a Lorentz factor less than 6,in agreement with its values for powerful 
Active Galactic Nuclei. 


Subject headings: Gosmology:theory,dark energy,SNe Ia,angular size-redshift re¬ 
lation,critical redshift,relativistic beaming 


INTRODUCTION 


There is now substantial observational evidence fjPeebles &: Ratrall2003l) that favors the 
existence of a smooth exotic cosmic component of energy of negative pressure. Going at 
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times under the name of a cosmological constant or quintessence or,at other times,dark 
energy,which we will adopt here,its true nature remains obscure. 

The unexpected faintness of high redshift type la supernovae (SNe I a)suggests that the 
universe is accelerating todavf Riess et ah 1998 : Perhnutter et ah 19991) .relentlessly driven 
by dark energy. When the SNe la results are combined with observations of the ampli¬ 
tudes of primordial fluctuations in the cosmic microwave background radiation the over¬ 
all picture seems to be one of an accelerating flat universe. Since the standard flat uni¬ 
verse, de spite its well-known s hortcomings,has long been favored on aesthetic and theoretical 
grounds flKolb fc Turneiill99Cll) .the hope has arisen that the injection of dark energy will cure 
its ills, particularly in regards of its age of the universe problem. Thus a major industry 
of investigating the constraints imposed by continuously updated astrophysical observations 
on the dark energy in rehned versions of the standard model has flourished in recent times. 
The present paper is one more contribution in this direction. 

When it is assumed that the dark energy,viewed in general to be time-dependent,does 
not interact with matter,the energy equations for nonrelativistic pressureless matter and 
dark energy decouple leading to conservation of matter and to the dark energy equation 
dpde/dz = 3(1 -1- pde/Pde)Pde, where Pde{z) and Pde{z) are the dark energy den¬ 

sity and pressure respectively and is the redshift. In this case a sol v able cosmologi- 


2002: 

Abdel-Rahman & Has 

rim 

2005)or a dehnite parameterization : 

or the equation o 

state 

w(z) = Pde!Pde is snggested ( 

Chevallier & Polarski 

2001: 

Bacciffahini 

2003 

Alam et ah 

2OO3I; 

Dicus & Renko 2 

004: Padmanabhan & Choudhurv 

2003 

Corasaniti et ah 

2004: Alam et ah 

2004; 

Johri 2004 

Johri & Rath 

200). 


Alternatively if one assumes that the dark energy interacts with matter (jOzeiill999l) .the 
energy equations for both are coupled and one needs a dehnite variation for pdeP'^ addition to 
specifying its equation of state. In this case the conservation of the energy-momentum tensor 
of the held equations,which holds when matter and the dark energy are noninteracting,is 
replaced by the conservation of the sum of this tensor and an extra appended tensorial 
piece representing the time-dependent dark energy. Here we follow this line: Specihcally we 
assume (a) pde ~ T,where T = p — 3p is the trace of the energy-momentum tensor of cosmic 
matter of density p and pressure p,and (b) a one parameter form for w{z). 


A variation A ^ T was introduced by iMaiernikI (120011.120031) for the cosmological con¬ 
stant A,the motivation being to identify the cosmological constant with a Lorentz-invariant 
scalar representing a form of quintessence. Thi s cosmology is reminiscen t o f similar earlier 
attem pts at identifying A with the Ricci scalar flAl-Rawaf fc Tahalll996alJbl : lAbdel-Rnhman 
19971) . The postulate A ~ T is interesting because it implies that the cosmological con- 






























































































3 


stant vanishes in the radiation-dominated cosmic era of flat cosmology so that the snccessfnl 
standard primordial nncleosynthesis predictions are nnaltered. In the matter-dominated 


era the postnlate rednces to A oc where H is Hnbble’s parameter. The cosmo 


constant variation A oc itself was widely discnssed in the literatnre flFreese et al 


Carva 


1997; 


Overdnin fc Cooperstock 


Vishwakarmal 1200 if) . In particniar 


ho et ah 1992 : Lima fc Carvalho 1994 : Arcnri fc Waea 1994 : Wetterichl 19951: 


ogical 


1987 


Arbab 


Carvalho et ah 


(jl992 l have pointed ont that it follows from dimensional argnments consistent with qnantnm 


gravity. Since snch argnments do not depend on the cosmological constant eqnation of state 
Pa = —Pa h is legitimate to regard them as eqnally valid for dark energy with w{z) ^ —1. 

Extending this postnlate to a dark energy with an eqnation of state of negative pressnre 
we take pde = where k is a dimensionless constant. A conseqnence of this is that the 
matter density parameter VLm is constant in the model. We take it to be 1/3. This is becanse 
a matter density parameter aronnd 0.30 seems to b e favored by ob servations indicat i ng that 
the dark energy acconnts for 2/3 of cosmic matter (jTnrnerll2002alJbl ). In fact iTnrnerl (I2002cr i 
has strongly argned a case for Vtm = 0.33T0.035 from measnrements of the physical properties 
of clnsters,CMB anisotropies and the power spectrnm of mass inhomogeneities. 


For the dark energy eqnation of state we consider 3 models with the one-parameter 
forms: {l)wde = w = const,{2)wde = —1 +wz,an.d {3)wde = —1 + w-r^, where w is constant. 


Model (1), viz, Wde = w = const < 0,is a generalization of the cosmological constant 
case Wde = —1- Strictly speaking a constant Wde is valid for the cosmological constant 
only. Yet models of cosmic evolntion driven by nonrelativistic matter and a qnintessence 


— 1), have been widely studied ( 

Ratra & Peebles 

1988; 

Chiba et al 

1997: 

Turner & White 

1997; 

Snereel & Pen 

1997; 

Frieman et al. 

198d; 

Caldwell et al. 

1998; 

Efstathiou 

1999 

Turner 

2002c 

1). In a number of these models (particularly those with tracking solutions),both the 


dark energy density parameter Qde 


flZlatev et al.lll999l:ISteinhardt et al. 


= SnGpde/SH"^) and Wde vary so slowly with redshift 


I 999 I: Efstathioul 1999 1 as t o instify the approximate use 
of ail effective equation of state parameter Weff ~ f Wang et al. 2000; Zhu et al. 


20041) . More generally,the absence of robust fundamental physics-based dark energy mod¬ 
els a nd the difficulty to observe a time de pendence of Wde from CMBRIjAnrich fc Steiner 
2003l )or from hts to luminosity distances (IDi Pietro fc ClaeskensI 120031) .admits the possi¬ 
bility of a Wde which is constant in some specihed range,and which arises as a model 


independen t approximation to the dark energy equation of state flKneller fc Strigaril 12003 


Cepal I2OO4J ) . (The cosmology with a dark energy ~ a ^ and de coupled from ordinary matter 


so that Wde = —1/3 has been rec ently discussed by one of us fj Abdel-Rahman! 120021) and by 
Abdel-Rahman fc Hashim ( 20051) ). 
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On the other hand parameterizations (2) and (3) are special cases of the two-parameter 
forms : Wde = Wr?p(0)+w^,and Wdp = W dp(0)+W j^ whi c h wer e proposed bv iHuterer fc Turner 


( 2001 ) and Weller fc Albrecht ( 20021) . and by Linder f 2003)respe c tively , and recently stud¬ 
ied,together with the case Wde = w = const, b MPicus &: Repkol (120041) . The form Wde = 
Wde{0) + wz diverges at very high redshifts whereas this difficulty is avoided i n the model 
Wdp = Wde{0) + Wj^, where Wde —^ Wde{0) + w a.s z ^ oo. But, as argued by iRiess et ah 


(120041) .a safer strategy, which we follow here,is to regard these parameterizations as only 
valid for low-z {z -C dec oupling redshift z ^^r) a nd describing the late behavior of dark en¬ 
ergy. This was done by iDicus &: Repkol (120041) who studied these parameterizations and 
found,on taking as prior = 0.3 in a flat universe,that a c onstant Wd^ = WdAO) = —1 is 


prefe r red by the fit to t he gold data for type la supernovae (IRiess et al.l |2004 iTonrv et ah 


20031: iBarris et ahl 120041) . Also it is already known that the cosmological constant scenario 


remains consist ent with tight con straints from new cosmic microwave background and galaxy 
clustering data (jMelchiorril 120041) . Quite genera lly observations seem t o require dark energy 
with present values Wde ~ —1 and Vide ~ 0.7 (jPeebles fc Ratral 120031) . With this in mind, 
and noting that recent SNela observations from HST do not indicate a rapid variation of 
Wde{z) away from its cosmological constant value,we pursue,for simplicity,the following ap¬ 
proach: we consider the preceding three Wde{z) parameterizations and set in them,a6 initio, 
Wde{^) = ~1- Then we investigate the constraints on the dark energy equation of state from 
recent supernova data and observations of the angular sizes of ultra compact radio sources. 

In section 2 we present the basic equations of the models. In sections 3 and 4 we examine 
the constraints on them from supernova and angular sizes data respec tively. In discussing 
angula r sizes we compare our results with those from a recent work by IJackson fc .Ta.nnetta 
(l2006l ).and also consider the implications of our angular size-redshift relations in the presence 
of relativistic beaming of the radio sources. Section 5 winds up the paper with a discussion 
of the results and some concluding remarks. 


2 . THE MODEL 

We consider a spatially flat FRW universe (a is the RW scale factor) 

= dt"^ — [dr^ -|- sin^ dd0^) (1) 

with cold matter of zero pressure and energy density pm and dark energy of density pde 
and pressure pde = WdePde- Denoting the scale factor today by oq,( subscript ”0” denotes 
present-day quantities),and dehning a/oo = (1 -|- z)”^,where z is the red-shift, Einstein’s 
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gravitational field equations can,in this case,be written as (a = S/SnG) 

a-^Pm + Pde) = H^ = HlE\z), 

3cr WdePm ‘2(^H 


(1 + ^Wde) 


= 


(1 + 3tCde)’ 


where q = —^ = 


z{l+z) 


2 is the deceleration parameter and H = a/a = 


(i+d 


( 2 ) 

(3) 

Hub¬ 


ble’s constant (an over-dot denotes time differentiation),with Hq = lOOh kms~^Mpc~^ = 
2.16hx being its present-day value {h is the normalized Hubble constant). Dehning 

the density parameters 

_ a ^pm ^ _ a ^Pde , .X 

we deduce from equation(l2]) that Qm + ^de = 1,valid at all times including t = tg. 


Combining equations 


|)we obtain 


(1+Z)dll^ 1 3Wde ^ X 1 3a ^Wdepde 


(5) 


In the Einstein-de Sitter(EdeS)standard model = 1 or pde = 0 so that g = 1/2. For 
the cosmological constant A case Wde = —1 so that q = |f2m —1 which admits an accelerating 
universe scenario provided flm <2/3. 

In this paper we assume that pde = i^T where T = p — 3v is the trace of the matter 
energy-momentum tensor and k a dimensionless constant flMajernikl 120011.120031 ). Then in 
the matter-dominated epoch of flat cosmolo gy we have fro r n equ ations (l2l) and (0]) that 
Pde = i^ Pm. = cind = const f)Maiernildl200ll. l2003l) . We further set,as was 

done bv iMaiernikI (1200 ll. l2003l) and argued in the introduction, = 1/3. We then obtain 
from equation(l5]), 

(1 -F z) dH^ 1 


q = -1 + 


2H^ 


dz =5 + '^'*- 


( 6 ) 


The rest of the paper investigates the consequences of this model for q using the different 
dark energy parameterizations discussed in the introduction. 


2.1. Parameterizations of Wde 

2.1.1. Model 1: Wde = w = const < 0. 
Inserting Wde = w = const < 0 in equation([6])yields 


E\z) = (1 + ;2)(3+2“’) 


(7) 
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2.1.2. Model 2: Wde = — 1 + wz., w = constant > 0 

Here equation(l6])shows that q > 0, {q < 0), for 2 ; > ^,{z < ^),implying a cosmic 
deceleration - acceleration transition at redshift = ^- In this case the solution of equation 
([6]) for if Ms ( 2 ; -C Zdec)- 

E\z) = (1 + exp {2wz). (8) 


2.1.3. Model 3: Wde = —1 + w 


(i+d’ 


w = const > 0 


The deceleration-acceleration cosmic transition occurs in this model at zt = 


1 


(2to-l) 


SO 


that we must have w > ^. In this case the solution of equation^ for if^ is( 2 ; -C Zdec)'- 


E\z) = (l + ;2)(^+2“’)exp [- 


2wz 


(1 + ^)- 


(9) 


3. TYPE la SUPERNOVAE 
3.1. The Distance Modulus 


For a flat universe the luminosity distance in units of Megaparsecs may be defined by 

dz 

dL = CH^\1 + z) ^ ^ c(l + z)M(0, ^), (10) 

Jo -nUJ 

where 

d{z„z,) = H,-\l + z,)-^ ( 11 ) 

E[z) 

In terms of di the predicted distance modulus is 

/rp = 51ogdL + 25. (12) 


We next obtain expressions for di and /ip in our models. In calc ulating /ip we use the 
widely accepted value for the Hubble constant Hq = 72kms~^Mpc~^ (IFreedman et al.ll2001 


Freedman &: Turneii 120031) . 


3.1.1. Model 1 


From equations(I71) and (ITOil 
di = c{l + z)‘^d{0, z) = 


2cHq ^(1 -|- ^) 


1 - (l + z)-’"-5 


1 + 2w 


w ^ -1/2, 
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di = c{l + zYd{0, z) = cHq^{ 1 + z)ln{l + z), w =—1/2. 
Hence by equation ([12]), 

2{l + z) 


(13) 


fip = 43.10+ 5 log I [l-(l + z)-"-2j|, w^-1/2, 

Up = 43.10 + 5 log [(1 + ;2) ln(l + 2;)] , tc = —1/2. 


(14) 


3.1.2. Model 2 


From equations (151) and flTUl 
di = c(l + zYd{D^ z) = exp(t(;)(l + z) 

where 


7 ( m; + ^,'«;(l + 2;) j -j (^w + ^,w 


'y{u,a) = / exp{—t)dt, Reu > 0, lar^fcrl < vr. 

Jo 


is the incomplete gamma function flAbramowitz &: StegnnI 119641) . Then 


so that 


exp{—wu)du 


-l+z 


/rp(2;) = 43.10 + 5 log I exp(t(;)(l + 2;) J 2 exp(— wm)^^ ) . (18) 


(15) 

(16) 

(17) 


3.1.3. Model 3 

Equal ions (Ej) and flTOl) give {w > 1/2) 

di = c{l + zYd{D.i z) = cHq^w^~^ exY){w){l + z) 

R 

r-l . 


1 W 


Hence 


= cHq ^ exp{w){l + z) / 2 Qxp{—wu)du. 


Pp(2;) = 43.10 + 5 log I exp(t(;)(l + 2;) J ^ exY>{—wu)du \ . 


(19) 


( 20 ) 













3.2. Supernova model predictions and observations 


3.2.1. Supernova observations 


Several astrophysics groups ljTonrv et al.ll2003l: iBarris et al. 


recent ly updated the original supernova data of iRiess et al 


199 


Riess et al. 

2004 

)have 

and 

Perlmutter et al. 


(119991) that provide d the hrst glimpse into an apparently accelerating universe. In particular 
Harris et ahl (j200J)have published photometric and spectroscopic observations of 23 super¬ 
novae in the redshift range 0.3396 < 2 ; < 1.031. Confronting our predictions for / ip with their 


data as analyzed by the BATM (Bayesian Adapted Template Match) method (jTonrv et al. 
2003l ).and calculated using Hq = 72kms~^Mpc~^, we have minimized with respect to the 
parameter w the statistic: 


x^(w) = f; 

i=l 


l^observed,i\ 


a: 


( 21 ) 


where the summation is over all 23 data points in Table 11 of iBarris et ahl (120041) and ai the 
corresponding uncertainties in the observed distance moduli. We discuss the application of 
this procedure to our 3 models. 


3.2.2. Model 1 


Using the first of equations (fTTD we calculated with the aid of equation (|2ID x^(tc). For 
w ^ —ll2 Figure 1 shows that the resulting curve has a minimum xlnin — 15.4 at ta = —0.7, 
with upper limits w = —0.41 and w = —0.28 at the 68% and 95% conhdence levels (c.l.) 
respectively (with 22 degrees of freedom- d.o.f.). For w = —1/2 corresponding to g = 0 
(coasting universe) we obtain,using the second equation in (fTT)) . xLn = 19-67. 


To discuss the implications of the value w = —0.7 for the age of the universe in this model 
we hrst note the following. For a hat universe with a Hubble constant Hq = 72kms~^Mpc~^ 
and contributions to the mass-energy density today of 1/3 and 2/3 of its total value from 
non-relativistic matter and dark energy respect ively,it is observed that the age of the universe 
is ISGyr with uncertainty of about ±1.5Gyr (iFreedman Sz Turnerl 120031) . A co nsistent age 


to = 14 ± O.bGyr is also determined from CMB anisotropy,independently of Hq (IKuox et al. 


20011 ). Moreov er,computer simulations o f Globular-cluster stars evolution produce ages of 


12.5 ± 1.5Gyr (iKrauss fc Chaboveiil200ll) . These estimates agree with values of to obtained 
by a variety of other m ethods, e.g. from r ates of cooling of old white dwarf stars or froin 
radioactive chronology (jOswald et al.lll996l) . Finally,assuming Wde = — 1 iTonrv et ahl (120031) 


deduce the constraint = 0.96±0.04,in agreement with the product (IFreedman fc Turner 
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2nn3h 

{Ho = 72 ± 8kms-^Mpc-^) x {to = 13 ± l.bGyr) = 0.96 ± 0.16. (22) 


The observed ages of the universe are therefore consistent with a consensus age of about 
13 ± 1.5Gyr ( Freedman fc Turner 2003 ). 


In the present model we have,from equation (171). 

r dz 2 

“ Jo {l + z)E{z) ~ 3 + ^- ^ 

Then w = —0.7 gives Hoto = 1.25. At the 68% c.l. w = —0.41 corresponding to Hoto = 
0.91, an estimate accommodated by equation ([22]). The coasting cosmology {w = —1/2) 
cor responds, as is well known,to Hoto = 1. 


3.2.3. Model 2 

In this model we have evaluated equation fl2T|) in the range 0 < w < 1 using equation flTSl) 
for the calculated distance modulus and plotted the results in Figure 2. We note that 
decreases monotonically as w increases from 0 but reaches a minimum xlnin — 16.5 at 
w = 1.1, corresponding to the transition redshift = ^ = 0.45. The value w = 2.55 
corresponds to the 68% c.l. limit. 


3 . 2 . 4 . Model 3 

Here we used equation ([201) in equation ([2Til and plotted versus w : 0.55 < tc < 2 in 
Figure 3. The curve has a minimum = 16.35 aX w = 1.7,corresponding to the transition 
redshift zt = 2 wZZi ~ 0.42. 

4. ANGULAR SIZE-REDSHIFT RELATION 
4.1. General formulae 

The angular size distance of a light source is 

dz 

dA{z) = d{0,z) = Hq\1 +z)-^ -^^ = {l + z)-^d{z), (24) 

where d{z) is the proper distance of the source. In a flat universe d{z) = aor{z), where r{z) 
is the source’s radial coordinate. 
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The angular size-redshift relation 6 = i/dA{z) where 6 is the source’s angular size,and i 
its intrinsic length,measured in j^arsecs {Ipc = 1.542xlO^^Geld“^)and assumed to be redshift- 
independent,is one of observational cosmology’s important tests of cosmological models. Like 
other classical kinematic tests it does not,generally,distinguish between cosmological models 
at low redshifts z -C 1 where the models are expected to converge. In fact for models with 
constant g,one has,for z -C 1, 


9 = 


iH. 


dA{z) 


l + -{3 + q)z + ... 


(25) 


which is formally the FRW result for small redshifts (jSandagelll988h . But for z > 1 there 
is less conhdence in the measurement s because of poss ible influences of poorly understood 
galactic evolutionary effects. However iKellermanI fll993l) has argued that ultra-compact radio 
sources with angular sizes (measured using VLBI: Very Long Baseline Interferometry) in the 
milliarcsecond {mas= 10“^ xl" = 4.8481x10“® radians) range (typica lly less than a hundred 
parse cs in extent) are deeply embedded in active galactic nuclei(AGNl fjPeterson 1997 : Krolik 
19981 1 and thus sheltered from extra-galactic evolutionary effects. Objects of this type have 
a fleeting existence (~ 100years),so it is reasonable to assume that characteristic parameters 
of their population (e.g. linear sizes)do not change on a cosmological time scale. 

Kellerman f 19931) showed that the angular size-redshift test for ul tra-compact sources 
favors the Einstein-deSitter flm = 1 canonical model. But subsequently iJackson &: Dodgson 


(119961.119971) demonstrated that the data is compatible with low-density constant-A mod- 

parameters for spat i ally fl at universes was 


els,indicating that the best choice of cosmologica 
Om = 0.2 and Oa = 0.8. In their latter work 


Jackson fc DodgsonI (119971) . utilized a data 


set of 337 ultra-com pact sources selected by GurvitsI ( 1994 ) from a 2.29 GHz survey by 


Preston et al.l ([l985|) comprising 917 sources with a correlated flux limit of approximately 0.1 
Jy (1 Jy = Jansky = Hz~^). From their study of this compilation they conclude 

that the canonica l mod el is ruled out by the observed angular diamete r-redshift rela tion. 
Later on Jackson (2004) rehned the analysis of Gurvits original data set ( GurvitsI 1994 ) and 
found for flat universes that Gm = 0.24 -|- 0.09/ — 0.07. Building on GurvitsI (Il994 ) ear¬ 
lier work iGurvits et ahl (jl999[) compiled a new data set of 330 compact radio sources which 
has,subsequently,been used by several a uthors in order to constrain the parameters of dif¬ 
feren t quintessence cosmological models (Vishwakarma 20^: TJma fc Alca^niz 2002; Alcani? 


2 OO 2 I: IZhu Sz Fuiimotoll2002l: iGhen fc Ratrall2003l: Ijain et aLll2003l) . 


We next give general formulae for the 6 — z relations in the present models. We write 
9 = where D = 6.87 x ih is the source’s characteristic angular scale (in mas). 
These expressions are (in mas): 
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Model 1: 


Model 2: 


D{w + i)(l + z) 

1 — (1 + ’ 

ln(l + z) ’ 


w ^ -1/2. 


w = -1/2. 


Z)w "'+2 exp(—w)(l + z) 

7 (w + |, w(l + z)) - 7 (w + w) 
D exp(—w)(l + z) 
exp(— 


(26) 

(27) 


(28) 


Model 3: 


9 = 


Dw^ 2 exp(—w)(l + z) 
D exp(—w)(l + z) 


f 1 2 exp(—wu)du 

1 + 2 

Using d'y{u,a)/da = a““^exp(—a) we find the small— 2 ; expansions of these equations: 
Model 1: 


(29) 


9 = — {1 + +3 + q)z + ... 


9 = -{l + ;-z+.. 


w ^ -1/2, 


w = -1/2, 


(30) 


(31) 


where from equation([6]) q = ^ + w. Thus both these equations formally coincide with 
eqnation fl25l) . 


Models 2 and 3: 


z 


1 + 7 + (8 + + - 


(32) 


where q = —^ + wz in model 2 and q = —^ + for model 3. In both cases the small- 2 ; 
expansion of 9 agrees with equation ([25]) on retaining only the ta-independent part of q. 


4.2. Critical redshift 

The existence of a critical redshift z^ corresponding to a minimum angular size can be 
qualitatively understood as follows. The reason for it is not just because,in the context of 
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cosmic expansion,was the light received today from a sonrce emitted when the sonrce was 
closer,but also,more importantly,because at a larger 2 ; the source of a standard linear size 
occupies a larger fraction of a large circle. From equations fl2Bl) - fl2^ the redshift Zm satisfies 
the equations: 


Model 1: 


Model 2: 


Model 3: 


Zm= W + 


- 1 , 


= e — 1 ~ 1.72, w = —1/2. 


(1 + = exp [w (1 + Zm)] 




(1 + = exp 


w 


1 + 


-1/2, 

(33) 

1 /2. 

(34) 

"2 exp{—wu)du. 

(35) 

3 

2 exp{—wu)du. 

(36) 


l + 2m 


Several authors f|Krauss fc Schramml Il993l: iLima &: A1 cani 71 l2000alJbl: IJain et ahl l2003l) 
studied critical redshifts in different models to find out how sensitive Zm is to variation of 
parameters like w. We address this question in §4.3.2. 


4.3. Constraints from angular size measurements 

analysis 

Our object now is to inve s tigate constraints on the parameters w and D using the data 
compilation of iGurvits et ahl fll999l) for the angular size measurements of milliarcsecond 
radio sources,observed at frequency z/ = 5GHz, in the redshift range 0.011 < 2 : < 4.72. 
The number of sources, originally 330, was reduced to those with a spectral index a in the 
range —0.38 < a < 0.18 and total luminosity L > 10‘^^W/Hz so as to minimize any possible 
dependence of angular size on a and also restrict the intrinsic size of the sources. These 
criteria were met by 145 sources which were then grouped in to 12 bins of 12-13 s ources per 
bin. This binned data was used in Figure 10 of the paper bv iGurvits et ahl fjl999l) . 


We attempted the determination of the best values of the models parameters D and w 
through a simultaneous minimization with respect to D and w of the statistic 
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X^(w,T>) = 

i=l 


measured 


a: 


(37) 
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where 9p^i denote the predicted angular sizes g i ven by equations ,Omeasured{zi) the 

corresponding observed values of iGurvits et ahl fjl999l) .and cxj the observation’s error of the 
sample’s ith bin. We have chosen the range of D to span the interval [0.10, 2.20]mas and w 
in intervals appropriate to each model. We adopt a minimization procedure,for 10 d.o.f,that 
produces minimum values for w values corresponding to different values of D in its 
range of variation. The least value among these minima,denoted XmmToi'responds to a 
minimization of x^ with respect to both D and w. The results are displayed in Tables 1, 2, 
and 3, and in Figure 5. The values of w and D corresponding to Xmin follows. 


Model 1 


Case w 7 ^ —1/2: 

Using equation fl2B]l in equation fl57D with w in the range [—1,0] we obtain (Table 1) 
Xmin = 4.72 corresponding to w = —0.45 and D = 1.25mas. Figure 4 shows, for —1 < ta < 0, 
the 68% and 95% conhdence contours in the w — D plane. 

Model 1 


Case w = —1/2 {Coasting cosmology): 

From equations fl27|) and flTTIl we deduce (see Figure 5 ) that Xmin = 4.75 at D = 1.3/q;29 
mas and D = 1.3lo;43 mas, where the errors produce the 68% and 95% conhdence limits 
respectively. 

Model 2 

Using equations (125]) and fl57|) and taking w in the range [0,3] we note,in this case,that 
X^ decreases monotonically from = 104.29 at D = O.lOmas and w = 3 to the minimum 
value Xmin = 4.47 at D = 1.45mas and w = 0.50 ( Table 2). The value w = 0.50 corresponds 
to zt = 1.00. The 68% and 95% conhdence contours in the w — D plane,for 0 < ta < 3, are 
displayed in Figure 6 . 

Model 3 

In this model equations (l2^ and fl37D reveal,for w in the range [0.55, 3], that descends 
from = 108.87 at D = O.lOmas and ta = 3 to a minimum Xmin = 4.58 at D = 1.40mas 
and w = 1.15 ( Table 3). This value of w corresponds to zt = 0.77. For 0.55 < w < 3, the 
68% and 95% conhdence contours are shown in Figure 7. 
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4.3.2. Calculation of the critical redshifts 


Equations fl33l) - fl36l) can be solved for Zm corresponding to specific values of w. The 
results are plotted in Figures 8-10 for models 1-3 respectively. The redshifts Zm for the best- 
£t w values are shown on the diagrams. We note that in both models 1 and 3 the b est-ht 
w value corresponds to Zm = 1.65. This reiterates the point made by Ijain et ahl (120031) 


that the minimum redshift test ” cannot by itself differentiate between different cosmological 
models” because different scenarios might correspond to the same Zm- The hrst and third 
models in this paper are generically quite distinct, yet they have the same Zm- In fact for all 
three models the values of Zm are quite close. In particular in model 2 the best-ht w gives 
Zm = 1-71 which almost coincides with Zm = 1-72 for coasting cosmology. Expectedly for 
re = 0 in model 1 we recover (Figure 8 ) the standard model result Zm = 5/4. 


4.3.3. 


Comparison with lCurvits et all U99ti) data 


The predictions of the 9 — z relations (I26l)- (I2^ for D and w corresponding to in 

each model are plotted,alongside the data of 

Gurvits et ah 

(199 

9),in 

Figure 11. ( Table 4 

gives log 2 ; and the corresponding log 6 * values for 

Gurvits et ah 

1999 

) data points and the 


end points log 6 ^- 1 - of their error bars). The inset magnifies the neighborhood of the minima 
of the curves. The curves,drawn for the best-ht values of w and cluster with in a narrow 


spread of log 9 and appear to be reasonably consistent with lGurvits et ahl (119991) data. None 


of the four models appears to be particularly favored over the others by this data. 


Lastly in the work of iGurvits et ahl (jl999l) a multi-parameter regression analysis of 
the data with ih a free parameter yields ^h = 23.8 ± 17.0pc. This corresponds to D = 
1.64 ± 1.17mas, a range within which fall all the D values in Table 5. 


4.4. Credibility of ultra compact radio sources as standard measuring rods 


In his 1993 analysis of compact radio sources iKellermanI (jl993l) found the angular sizes 
to be essentia l ly ind ependent of redhift in the interval 0.5 < z < 3. In the latter work of 
Gurvits et ahl (jl9991 1this featur e appeared to per sist for median angular sizes of sources with 
z > 0.5. Although according to IKellermanI (jl993l) .and as mentioned earlier on,it is reasonable 
to assume that milliarcsecond ultr a-compact rad i o sour ces are not ahected by evolutionary 
ehects,it is still important to treat iGurvits et ahl (119991) data ,as emphasized by them, with 
caution. A widely-accepted model of an ultra compact radio source is one in which a central 
low luminosity object is straddled by a pair of radio-bright lobes emitting synchrotron radia- 
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tion. The lobes are sustained by two hot gas iets f|Peacocklll999l: |Petersonlll997l: lKrohklll998l) . 
In such a picture differences in the spectral index between the core,the central engine, and 
its jet com ponents may intr oduce a ’’linear size-redshift ” dependence even in the absence of 
evolution. iFrev et ahl fll997l) have however argued that such a dependence is likely to be weak. 
Nevertheless more data at various frequencies going beyond limited source samples can,in 
particular,help verify the importance of such a dependence and also,more generally,enhance 
further the c redibility of milliarcs e cond radio-sources as standard measuring rods. With this 
goal in mind I Jackson fc .Tannettal fl200^have recently studied the implications of an updated 
sample of sources in the IPreston et ah f 1985 ')catalogue. We shall shortly discuss t heir work 
and compare their results to ours. Another aspect of interest is relativistic bearning (jPeacpck 


19991: lPetersorull997l : lKrolikill998l) .the process by which the Doppler effect modihes the ap¬ 


pearance of a radio source with lobes. For a lobe whose jet axis is oriented close to the line of 
sight from the source to the earth relativistic beaming generally increases the apparent radio 
power of the source and de creases its ang ular size. The implica tions of this orientation bias 
on the FRW 9 — z relation fjSandagelll988l )have been studied bv lDabrowski et al.l (119951) with 
the aim of evaluating the statistical confidence in the observational data,p articularly in rela¬ 
tion to estimates of the value of go- A similar study was also undertaken bylStepanas fc Saha 
(119951 ). We shall consider a simple variation of the iDabrowski et ahl (ll995|) approach in order 
to assess the impact of this orientation bias on our calculations. 


4.4-1- Comparison with the work o uJackson & Jannettd 12009) 


Jackson fc .lannettal ([200^ reexamined the lPreston et al.l (1 19851 ) catalogue of ultra com¬ 
pact radio sources,updating it with respect t o both r e d shi ft and radio information and 
replacing the original choice of 337 sources by iGurvitsI (119941) with a sample of 613 objects 
in the red shift range 2 ; = 0.0035 to 2 ; = 3.787. 


Over its red shift range of 0.003 < 2 ; < 3.8 iGurvitsI (119941) 337 objects sample contains 
sourc es of luminosit y that varies over three orders of magnitude,from 0.01 to 1 Jy (see Fig 
3 in (IGurvitsI 119941) j. Moreover the sources in the high red shift range ^ > 0.5 exhibit a 
smaller dispersion in luminosity as compared to sources with 2 : <0.5. T h is sug gests that 
the higher red shift object s have similar linear sizes. This is why IGurvitsI (jl994j) . and also 
Jackson fc Jannettal (120061 ) .ignore objects with 2 : < 0.5. 


This mismatch betwe en the luminos i ty of h igher and lo wer z sources is also a pparent in 
the 330 sources sample of iGurvits et al.l (119991) (see Fig 3 in iGurvits et al.l (ll999l D. This is 
why that we also limit ourselv es in this paper to lGurvits et al.l (Il999l) data with > 0.5 (see 
Fig 10 in iGurvits et ahl (jl999l) and Table 4 in the present paper). 
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Of the 613 sources in the sample of l.Tacksnu fc Januettal fl2006h 468 have 2 : > 0.5. These 
were placed in 6 bins with 78 objects in each bin. But IJackson fc .lannettal (120061) show that 
the cosmological parameters which best £t the 6 — z data in this case,including a best-fitting 
value i = 7.7bh~^kpc for the source’s intrinsic length,are not very se nsitive to the choice 
of bin ning. Rather than relying too heavily on a particular binning IJackson fc .la.nnetta 
(120061) produce data points corresponding to bin sizes of 76, 77 and 78 objects and take 
for their {z, 9) points values that are a composite of these three cases. Their resulting data 
points are: (z, e{mas)) = (0.6153,1.4624), (0.8580,1.2801), (1.1527,1.1599), (1.4200,1.1448), 
(1.8288,1.1760), and (2.5923,1.2374),with standard de viation a = 0.00603 . Us ing this 
data,and marginalizing over the sources intrinsic length Jjackson fc .lannettal ( 120061 ) obtain 
the best-£t parameter values VLm = 0.25-1-0.04/ —0.03, IIa = 0.97-1-0.09/ —0.13 and the curva¬ 
ture parameter k = 0.22-1-0.07/ —0.10,where the error bars are 68 % con fidence limits. These 
result s do not agree with our models. In fact the log 6 — log z curve fo rl-Tackson fc .Tannetta 
(120061) results would lie well below the base of Figure 11. We have con firmed this disagree¬ 
ment by doing a analysis of the predictions of our models against IJackson fc .Ia.nnetta 
(l2006l) angular size results. We will comment on this matter in the summary and conclusions 
section. 


4-4-2- Effect of orientation of lobe jet axis on qo:Method of Dabrowski et al. (199^) 


Dabrowski et ahl (119951 ) gave qualitative and quantitative discussions of the effect of 
relativistic beaming on the FRW 6 — z relation. They consider a simulation sample of 
compact radio sources,with each source composed of two identical but oppositely directed 
jets having an angle 0 to the line of sight,where 0 < 0 < vr. Each source emits an isotropic 
power-law spectrum with flux density 


s = ± 

di 


3+0 


7(1 — f3 cos(/)) J ’ ^ ^ 

where L is the source’s luminosity, di is the luminosity distance, a is the spectral index, is 
the jet speed/c,and 7 = (1 — is the Lorentz factor of the jets. For a quasi-continuous 

j et formed out of finite-lifetime blobs the appropriate power index above is reduced to ( 2 -|-a) 
( Peacock 19991) . 


The inclusion of the source’s lobe orientation in the FRW 9 

m,ql{l + zf 


z relation 


qaz + (90 - 1 ) (( 2*9 + 1 )“'" - 1 ) 

replaces the source’s intrinsic length f sin 0 where is the projection of f in the 

plane of the sky (£±(0 = f) = ^ source whose jet axis lies along the plane of the sky) 
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It follows from equ ation fl38|l that the maximum allowed angle in the sample is given by 
flDabrowski et al.lll995l) : 


(j)max{z) = arccos 
where Sum denotes the sources flux density limit. 



Slimd\ 


-1 

3 +q 


Dabrowski et ahl (jl995h dehne the average lobe orientation angle by 


4>{z) = 


r 0 max • 111 

Jq sm (paq) 


(40) 


(41) 


Since 0 is small on e can assume i i jz ) = and insert (-i_{z) into the 6 — z relations. 

The essence of the Dabrowski et ah f 19951) quantitative method is to investigate the effect 
of this modihcation on the predictions of equation 


4 . 13 . 


Optimized relativistic beaming: an alternative approach to \Dahrowski et al\ h99d) 

calculation 


The twin-jet model of iDabrowski et ahl (119951 ) assumes that the measured angular size 
of the source c orresponds to the separation of the jets projected onto the plane of the sky. 
JacksonI (120041) considers this feature to be unrealistic because the flux of the receding (away 
from us) jet is very much smaller than that of the advancing(towards us)jet,for example 
by a factor of up to 10^ for a jet Lorentz factor 7 of 5, and may therefore be neglected 
( Lind fc Blandford 1985 ). 


The relativistic beaming hypothesis is primarily based on two parameters: the Lorentz 
factor 7 (or equivalently the parameter = the jet velocity/c ) and the viewing angle 0. In 
a study of 25 quasars with 2 ; between 0.1 and 2.5 (some quasars with their nonstellar radi¬ 
ation spectrum and evidence of jets and ext ended emission features also display relativistic 
beaming effects iVermeulen fc CohenI (jl994l) hnd that 7 in creases with 2 : if qn = 0.05 but 
remains constant at 7 ~ 10 if go = 0.5. With this in mind iDabrowski et ahl (119951) take, in 
their analysis of a simulation data, 7 = 10 , although they point out that this choice does 
not qualitatively affect their results. 


However a choice of 7 linked to a positive go,as in the 1995 paper bv iDabrowski et ah 


(119951) .is hard to justify in these days where go < 0 (accelerating universe) is very probable. 
Moreover the re are other reasons that might render an arbitrary choice of 7 fraught with 
uncertainties (IUbachukwu fc Chukwudel l 2002 l) . In addition, geometric projection effects of 
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components of jets of galaxies or quasars can lead to dramatically misleading values of 7 . 
For example the transverse motion of blobs of gas across the sky in a well-studied jet in 
the galaxy M87 ( at about ISMpc away from us ) appears to be superluminal corresponding 
to a velocity nearly 6 times the speed of light! fIChaisson fc McMillanI l2002l) . This cosmic 
illusion of faster-than-light velocities is simply a projection effect, pr oduced by the blobs of 
gas m oving at a very small angle 0 at a near-velocity-of-light speed (jChaisson fc McMillan 
2 OO 2 I). It is worthy of mention that Doppler boosting was hrst discussed by IShklovskv 


f 1964a.bi yjto explain the apparently one-side jet inM87. 


In this paper we choo se to relate 7 to 0 by opti r nizing the re l ativis tic b eaming effect. _ 

Iris t echnique was used bv IVermeulen fc Cohen! fl 19941) iUbachukwul fjl999l) .and lUbachukwu fc Chukwude 


fl2nn2[i 


Dehne the Doppler factor which fundamentally characterizes relativistic beaming by 
(see equation (l38|) j: 

5 = [ 7(1 — 0cos0)]“^ (42) 


Optirnization of relativistic beaming is obtained by requi ring H = 0. This yields f Vermeulen &: Cohen 
1994 IUbachukwulll999l: lUbachukwu fc Chukwudell 2002 l) 0 = arccos 0, or 0 = arcsin Since 

the lobe orientation bias replaces i by f'sin0 in the 6 — z relation, D = 6.87 x in 

equations 


is replaced by —. In what follows we discuss the constraints of the obser¬ 
vational data on —. 

7 

Clearl y application o f the y ^procedure of §(4.3.1.) to the modihed 6 — z relations using 
the binned iGurvits et ahl fll999l) data assigns now the best-£t values of D to ^ and leaves 
those of w (and hence of q'o)undisturbed. Taken together the three models imply (see Table 5) 


1.25 < — < 1.45 mas, or, equivalently, 18 h “^7 < i < 21 h “^7 pc. Since milliarcsecond ultra¬ 
compact radio sources are typically less than a hundred parsec in extent the last inequality 
would imply that 18h ~S < i < lOOpc, or 7 < < 6 . (If we use h = 0.72 we get 


7 < 4). It is known (lAbramowicz et al.lll990l: iPeacockI Il999l )that the most powerful AGN 
jets seem to have Lorentz factors of the order of these upper bounds. 


5. SUMMARY AND CONCLUSIONS 

We have studied in this paper supernova and compact radio sources angular sizes con¬ 
straints on three cosmological models whose dynamics is driven by non-relativistic matter 
of density parameter VLm = | and a smooth time-dependent dark energy component with 
density pde oc T = p — 3p (T is the matter energy-momentum tensor) and equation of 
state Wde = being either constant or redshift-dependent, with Wde{z = 0 ) = —1 
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in the latter case. The variation pde ~ T implies that the dark energy vanishes in the 
early nniverse leaving the standard model’s primordial nucleosynthesis predictions intact. 
This is somewhat reminiscent of the modihed general relativity model of lAl-Rawaf fc Taha 
fll996al]bl) which can be cast in the form of a variable-A cosmology with A oc R, R being 
the Ricci tensor. There the radiation density pr ~ a~^ ~ in the early universe, as in 
standard flat cosmology. However the Friedmann equation is modihed so that the standard 
cosmic expansion rate is altered by the factor |(1 + i^/rf)) where rj, a constant, satishes 
0<ri<l,ri = lm general relativity. Here the postulated dark energy ansatz does not 
affect the standard cosmic expansion rate of the early universe. It reduces in the matter- 
dominate d phase to php ^ H^,a. variation that was extensively studied for the cosmological 


constant f Overdnin fc Cooperstock 1998h . Our main results are summarized in Table 5. 


For the hrst mode l with Wdp constant (implying by equation ([ 6 ]) that q = qo = constant) 
the supernova data of iBarris et al.l (120041) admits an accelerating universe with Wde = —o.7 
and go = — 0.2 = 15.4) and a coasting universe with Wde = — 19.67). On 

the other hand iGurvits et ahl fjl999l) angular sizes data allows a mildly decelerating universe 
with Wde = —0.45,go = 0.05 and characteristic angular scale D = 1.25mas {Xmin = 4.72) 
and also a coasting cosmology with Wde = —0.50 and D = l.SOmas {Xmin — 4.75). The 
value D = l.Smas corresponding to a source’s intrinsic length i = 18.92h~^p c is ver y close 
to D = 1.28 obtained as the best-£t value in an earlier model by I Jain et al.l (120031) where 
the scale factor is essentially linear in t,viz^a ~ 

The results of the constant-w^e model (and also those of the other models in this paper) 
do not extrapolate to the radiation dominated universe since there,as we have pointed out in 
the introduction, pde vanishes. Nevertheless an accelerated expansion in model 1 may pose 
problems for structure condensation after matter dominance: structure tends not to form in 
the presence of cosmic acceleration. However the modest acceleration (g = go = —0.2) in 
model 1 corresponds to Wde = —0.7 with upper limits Wde = —0.41 and Wde = —0.28 at the 
68 % and 95% conhdence levels respectively. These values of Wde correspond from equation 
([ 6 ]) to g = 0.09 and g = 0.22 respectively so that within the limits of th e quoted confidence 
levels a decelerating universe is not excluded. In fact it has been argued (jVishwakarmall2003l ) 
that absorption by intergalactic dust of light travelling over immensely long distances might 
explain the faintness of extragalactic SNe la obviating the need for a cosmic acceleration 
based explanation. Yet we are inclined to believe that our supernova results for t he coasting 


universe scenario are more robust because of their concordance with those of Jain et ah 


(120031) which are obtained by a different approach. A flat cosmology with dark energy of 
constant Wde and ma tter density parameter held at Qm = 0.3 was also considered by 
Dicus &: Repkol (1200411 who however assumed noninteracting dark energy and used the larger 
set of the supernova ’’gold” data. They hnd in this case the preferred (cosmological) constant 
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value Wde = — 1- 


Models 2 and 3,which by construction describe a present-day accelerating universe with 
go = —0.5,are consistent with the data for universes that undergo deceleration-acceleration 
transitions at redshifts in the range 0.42 < < 1- In the standard ACDM cosmological 

model zt = 0.67 for Q m ^ 0.30, in c ontra. st with the observational value zt = 0.46±0.13 from 
the SNela analysis of iRiess et al.l (120041 ). Predictions of zt in seven popular qui ntessence 
model s inspired by supergravity or Ml string theory have recently been studied bv [Gardner 
(I 2 OO 5 I) who noted that all of them can, in the low-z approximation Wde ~ Wo + WiZ,{Q < z < 
5),mimic the ACDM model. Here in models 2 and 3 the SNela data give Z;: = 0.45 and zt = 
0.42,both values being very close to the observational result zt = 0.46 o fiRiess et al.l fl2004r) . 


On the other hand consistency of the models with the angular size data flGurvits et ah 


199911 


yield zt = and zt = 0.77, to be compared with zt = 0.67 from the ACDM. Inspection 
of Figure 11 and Table 5 lead s to an import a nt con clusion:The model’s angular size-redshift 
curves drawn in Figure 11 fit iGurvits et ahl (119991 1 data equally well. This is a reflection of 
the small values and span of Xmm seen in Table 5 for the results : for the four models 
Amin varies only from 4.47 to 4.75. By comparison in the supernova test the values 
are considerably larger and he in the range 15.4 to 19.67 so that one has more faith in the 
agreement of the models with the angular size data. 

The predictions of the models for the critical minimum redshift in the angular size- 
redshift relation give Zm values in the narrow interval [1.65,1.72] compared to the standard 
model result Zm = 1.25. The constant-tCrfe 7^ ~l/2 and the third models have the same z^ = 


1.65 w hereas the coasting universe and model 2 have ZmX.72 (the same as that of I Jain et ah 
(j2003l) Iand Zm = 1-71 respectively. Thus the minimal redshift cannot,by itself,effectively 
discriminate between these models. 

We have also compared the predictions of our mode ls with those o f the recent work 
of Jackson fc Jannetta (2006)in which they revisit the old Preston et ah (jlOSh lcatalogue of 
ultra compact radio sources and reconsider an angular size-redshift data set in the light 
of modern preferences of the cosmological parameters. We hnd that their results do not 
agree with our models. Ii i a sense this is not surpri sing since the underlying premises of 
the two works are different ! Jackson fc Jannettal (120061) approach is based on using one simple 
potential to test the hypothesis that vacuum energy is constant. The present models are 
phenomenological and based on a time-depend ent dark energ y coupled to matter. Hence 
there is no overlap between the two approaches (I JacksonI 120071) . 


Finally,we tested our angular size-redshift relations in the presence of relativistic beam¬ 
i ng. R elativistic beaming in the FRW Q — z relation was investigated by iDabrowski et ah 


(119951) . Here we used a simple variation of their method in a two-jet model of radio sources 
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with the advancing (towards uslie t axis close to the line of sight. In this pictnre consistency 
of the models with iGurvits et ahl fll999l) data is fonnd for a beaming Lorentz factor 7 < 6 . 
Snch an npper bonnd is consistent with valnes of the Lorentz factor for powerfnl AGN jets 
flAbramowicz et ahlll990l: lPeacocldll999l) . 


To conclnde, we have presented dark ene rgy models that are in reasona ble agreement 


with t he snpernova data of lBarris et ahl (120041) and in good agreement with the iGnrvits et ah 


(119991) compact radio sonrce angnlar size versns redshift binned data. The three mode l s that 
we have stndied are simplified versions of ones recently considered by iDicns fc Repkol (120041) 
in a different context. Bnt as remarked by these anthors comparing models for the eqnation 
of state of dark energy will remain something of a mng’s game nntil there exists snbstantially 
more data at higher redshifts. 
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Table 1. — w — D values for Model 1 


D 

w 


D 

w 


0.1 

0 

106 

■ 

1.25 

-0.45 

4.72 

0.2 

0 

85.44 

■ 

1.30 

-0.55 

4.74 

0.3 

0 

66.53 

■ 

1.35 

-0.60 

4.75 

0.4 

0 

50.16 

■ 

1.40 

-0.65 

4.80 

0.5 

0 

36.33 

■ 

1.45 

-0.70 

4.88 

0.6 

0 

25.03 

■ 

1.50 

-0.80 

4.97 

0.7 

0 

16.27 

■ 

1.55 

-0.85 

5.08 

0.8 

0 

10.04 

■ 

1.60 

-0.90 

5.20 

0.9 

0 

6.36 

■ 

1.65 

-0.95 

5.35 

0.95 

0 

5.47 

■ 

1.70 

-1 

5.51 

1.00 

-0.05 

5.17 

■ 

1.80 

-1 

6.28 

1.05 

-0.15 

5.02 

■ 

1.90 

-1 

7.96 

1.10 

-0.20 

4.90 

■ 

2.00 

-1 

10.54 

1.15 

-0.30 

4.80 

■ 

2.10 

-1 

14.02 

1.20 

-0.40 

4.76 

■ 

2.20 

-1 

18.41 
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Table 2. — w — D values for Model 2 


D 

w 


D 

w 


0.1 

3 

104.29 

■ 

1.25 

1.00 

5.00 

0.2 

3 

80.97 


1.30 

0.90 

4.78 

0.3 

3 

60.91 


1.35 

0.70 

4.61 

0.4 

3 

44.12 


1.40 

0.60 

4.50 

0.5 

3 

30.59 


1.45 

0.50 

4.47 

0.6 

3 

20.33 


1.50 

0.40 

4.52 

0.7 

3 

13.32 


1.55 

0.30 

4.64 

0.8 

3 

9.58 


1.60 

0.20 

4.84 

0.85 

2.80 

8.85 


1.65 

0.15 

5.07 

0.9 

2.50 

8.24 


1.70 

0.10 

5.39 

0.95 

2.20 

7.65 


1.75 

0 

5.78 

1.00 

1.90 

7.08 


1.80 

0 

6.28 

1.05 

1.70 

6.56 


1.90 

0 

7.96 

1.10 

1.50 

6.09 


2.00 

0 

10.54 

1.15 

1.30 

5.66 


2.10 

0 

14.02 

1.20 

1.10 

5.30 

1 

2.20 

0 

18.41 









Table 3. ~ w — D values for Model 3 


D 

w 


D 

w 


0.1 

3 

108.87 


1.25 

1.85 

4.83 

0.2 

3 

89.00 


1.30 

1.55 

4.70 

0.3 

3 

71.25 


1.35 

1.35 

4.61 

0.4 

3 

55.63 


1.40 

1.15 

4.58 

0.5 

3 

42.13 


1.45 

0.95 

4.61 

0.6 

3 

30.77 


1.50 

0.75 

4.69 

0.7 

3 

21.53 


1.55 

0.65 

4.82 

0.8 

3 

14.42 


1.60 

0.55 

5.04 

0.9 

3 

9.43 


1.70 

0.55 

6.28 

1.00 

3 

6.57 


1.80 

0.55 

8.59 

1.05 

3 

5.94 


1.90 

0.55 
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1.10 

2.65 
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1.15 

2.35 

5.28 


2.10 

0.55 

21.95 

1.20 

2.05 

5.03 


2.20 

0.55 

28.55 


Table 4. 


Angular size-redshift data fiGurvits et al.l 1199911 



logz 

logO 

a 

loge+ 

1 

O 

1 

-0.29 

0.785 

3.216 

0.961 

0.461 

2 

-0.20 

0.703 

2.346 

0.840 

0.420 

3 

-0.12 

0.585 

1.504 

0.722 

0.374 

4 

-0.05 

0.550 

3.170 

0.823 

-0.427 

5 

0.04 

0.597 

1.277 

0.704 

0.411 

6 

0.10 

0.550 

2.683 

0.765 

-0.075 

7 

0.15 

0.562 

0.634 

0.621 

0.482 

8 

0.19 

0.332 

1.354 

0.537 

-0.107 

9 

0.29 

0.389 

0.895 

0.508 

0.196 

10 

0.34 

0.690 

1.078 

0.768 

0.573 

11 

0.41 

0.525 

0.734 

0.623 

0.433 

12 

0.55 

0.633 

0.990 

0.731 

0.535 
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Table 5. Summary of the results 


Model 


SNela 



Angular size 


Critical redshift 


w 

A.min 

Zt 

Qo 

w 

D 

(mas) 

Y^ ■ 

A. mm 

Zt 

Qo 

^m 

1 

-0.7 

15.4 

_ 

-0.2 


1.25 

4.72 

_ 

0.05 

1.65 


-0.5 

19.67 

- 

0 


1.30 

4.75 

- 

0 

1.72 

2 

1.1 

16.5 

0.45 

-0.5 


1.45 

4.47 

1 

-0.50 

1.71 

3 

1.7 

16.35 

0.42 

-0.5 

1.15 

1.40 

4.58 

0.77 

-0.50 

1.65 
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Fig. 1.— Model 1 (g 7 ^ 0): Plot of equation versus w. 
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Fig. 2. — Model 2: Plot of equation fl2l|) .Y^ versus w. 
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Fig. 3.— Model 3:Plot of equation fl^ versus w. 
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Fig. 4.— Model 1 (g 7 ^ 0) : Confidence contours in the w — D plane. 
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Fig. 5.— Model 1 (g = 0): Plot of equation (ETj), versus D. 




























40 



Fig. 

bars 


11.— Predicted 6 — z cur ves for the models. Also shown are the data points and error 
from iGnrvits et ahl (119991 ). 

































